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Abstract 

We present some observations on a restricted variant of unitary Cayley graphs 
modulo n, and implications for a decomposition of elements of symplectic operators 
over the integers modulo n. We define quadratic unitary Cayley graphs Gn, whose 
vertex set is the ring ^72, and where residues fi, b modulo n are adjacent if and only 
if their difference is a quadratic residue. By bounding the diameter of such graphs, 
we show an upper bound on the number of elementary operations (symplectic scalar 
multiplications, symplectic row swaps, and row additions or subtractions) required 
to decompose a symplectic matrix over • also characterize the conditions on 
n for Gn to be a perfect graph. 

1 Introduction 

For an integer n ^ 1, we denote the ring of integers modulo n by Z„, and the group of 
multiplicative units modulo by Z^. A well-studied family of graphs are the unitary 
Cayley graphs on Z^, which are defined by X„ = Cay(Z„,Z^). These form the basis of 
the subject of graph representations [1], and are also studied as objects of independent 
interest: see for example [2-5]. 

We consider a subgraph Gn ^ A'„ of the unitary Cayley graphs, defined as follows. 
Let Qn = {u^ I u G Z^ } be the group of quadratic units modulo n (quadratic residues 
which are also multiplicative units), and T„ = ±Qn. We then define Gn = Cay(Z„,T„), 
in which two vertices given by a, 6 G Z„ are adjacent if and only if their difference is a 
quadratic unit in Z„, i.e. if a — 6 G |±ti^ | m G Z^}. In the case where n = 1 (mod 4) 
and is prime, G„ coincides with the Paley graph on n vertices: thus the graphs G„ are a 
circulant generalization of these graphs for arbitrary n. We refer to Gn as the (undirected) 
quadratic unitary Cayley graph on Z„. 



* niel. debeaudrap@gmail. com 



We present some structural properties of quadratic unitary Cayley graphs G^. In 

particular, we characterize its decompositions into tensor products over relatively prime 
factors of n, and categorize the graphs Gn in terms of their diameters. From these results, 
we obtain a corollary regarding the decomposition of symplectic matrices S G Sp2^(Z„) 
in terms of symplectic row-operations, consisting of symplectic scalar multiplications, 
symplectic row-swaps, and symplectic row- additions/subtractions. We also characterize 
the conditions under which quadratic unitary graphs are perfect, by examining special 
cases of quadratic unitary graphs which are self-complementary. 

Notation. Throughout the following, n = p^^p^^ ■ ■ -p^* is a decomposition of n into 
powers of distinct primes, and o" : Z„ — )> Zp^mj © ■ ■ ■ © Zp^m^ is the isomorphism of rings 
which is induced by the Chinese Remainder theorem. (We refer to similar isomorphisms 
p : Z„ — > Zm © Ztv for coprime M and N as natural isomorphisms.) We sometimes 
describe the properties of Gn in terms of the directed Cayley graph = Cay(Z„,(5n), 
whose arcs a ^ b correspond to addition (but not subtraction) of a quadratic unit to a 
modulus a e Z„; we may refer to this as the directed quadratic unitary Cayley graph. 

2 Tensor product structure 

By the isomorphism Z^ = Z^^mi © • • • © Z^^mt induced by a, unitary Cayley graphs X„ 
may be decomposed as tensor products Xn = Xp^n^ 0- • •(S)Xp^mt of smaller unitary Cayley 
graphs (also called direct products [5] or Kronecker products [6], among other terms): 

Definition I. The tensor product A^B of two (di-)graphs A and B is the (di-)graph with 
vertex-set V{A) xV{B), where {{u,u'), {v,v')) G E{A(S)B) if and only if {{u,v), {u',v')) e 
E{A) X E{B)} 

Corollary 3.3 of [5] gives an explicit proof that Xn — Xp^m^ ® • • • ® ^pt™*; a similar 
approach may be used to decompose any (di-)graph Cay(i?, M) for rings R — Ri®- ■ -^Rt 
and multiphcative monoids Mi — Mi © • • • © where Mj C Rj. For instance, as 
Q„ = Qpi""! © ■ ■ ■ © Qpt"^t , it follows that r„ = Fp^mi © ■ ■ ■ © Fp^mt as well. 

It is reasonable to suppose that the graphs Gn will also exhibit tensor product struc- 
ture; however, they do not always decompose over the prime power factors of n as do 
Xn and F„. This is because T„ may fail to decompose as a direct product of groups over 
the prime-power factors p^\ By definition, for each j, we either have Tp.mj — Qpruj or 
Tp^.^j = QpP^i © (—1); when Qp^.^j < ^p^^i for multiple pj, one cannot decompose r„ over 
the prime-power factors of n. We may generalize this observation as follows: 

Theorem 1. For coprime integers M,N 1, we have Gm (8) Gat — Gmn if o-nd only if 

either -1 G Q\i or -1 G Qy. 

^We write Ai (g) {A2 (g) A3) = {Ai (g) A2) X .4,3 — Ai A2 ® A3, and so on for higher-order tensor 
products, similarly to the convention for Cartesian products of sets. 



2 



Proof. We have Gm 8) = Gmn if and only if Tm ® TV = Tmn- Let p : Zmat — > 
Z^vf ©Z^v be the natural isomorphism: this induces an isomorphism Qmn — Qm®Qn, and 
will also induce an isomorphism Tmn — Tm®T]s[ if the two groups are indeed isomorphic. 
Clearly, a{TMN) ^ Tm ® T^; we consider the opposite inclusion. 

If —1 ^ Qm and —1 ^ Qat, we have (—1, 1), (1, —1) ^ Qm © Qiv; as both tuples are 
elements of Tm ® T^v, but neither of them are elements of ±{Qm ® Qn) — (^{'^Qmn) — 
<j{Tmn), it follows that Tmn and Tm © TV are not isomorphic in this case. Conversely, 
consider -u G Z^ arbitrary, and let {um,un) = p('w). If —1 G Qm, let i e Zm such that 
= —1: for any sm, sat G {0, 1}, we then have 



Remark. The above result is similar to [8, Theorem 8], which uses a "partial transpose" 
criterion to indicate when a graph may be regarded as a symmetric difference of tensor 
products of graphs on M and N vertices; the presence of —1 in either Qm ot Qn is equiv- 
alent to Gmn being invariant under partial transposes (w.r.t. to the tensor decomposition 
induced by p). 

Corollary 1-1. For 1, let n — p^^ ■ ■ ■p^'^N he a factorization of n such that pj = 1 
(mod 4) for each 1 ^ j ^ t, and N has no such prime factors. Then Gn = Gp^m-^ (g) • • • (g) 

Gp^m-r (g) Gn ■ 

Proof. For pj odd, Z^mj is a cyclic group [7] of order (pj — l)p^^ ^ in which —1 is the 
unique element of order two: then —1 is a quadratic residue modulo p^^ if and only if 
Pj = 1 (mod 4). As this holds for all 1 ^ j ^ r, repeated application of Theorem 1 yields 
the decomposition above. □ 

Corollary 1-2. For 1, we have Gn — Gp-^m^ (g) • • • (g) Gp^mt if and only if either n has 
at most one prime factor pj ^ 1 (mod 4), or n has two such factors and n = 2 (mod 4). 

Proof. Suppose that Gn decomposes as above. Let N be the largest factor of n which 
does not have prime factors p = 1 (mod 4): we continue from the proof of Corollary 1- 
1. By Theorem 1, Gn itself decomposes as a tensor factor over its prime power factors 
p^^i^ , . . . ,pp if and only if there is at most one such prime pj such that — 1 ^ Qp/^j- 
However, by construction, all odd prime factors pj of N satisfy pj = 3 (mod 4), in which 
case —1 ^ Qp/^j for any of them. Furthermore, for m ^ 2, we have r e Q2m only if 
r = 1 (mod 4); then -1 E Q2m if and only if 2"" = 2. Thus, if G„ = Gp^m, © • • • (g) 
Gp^mt, it follows either that = p"^ for some prime p = 3 (mod 4), in which case 
the decomposition of Corollary 1-1 is the desired decomposition, or N = 2p"^ for some 
prime p = 3 (mod 4), in which case n = 2 (mod 4). The converse follows easily from 
Corollary 1-1 and Theorem 1. □ 





(1) 



Thus Tm ® Tjv ^ cr{TMN); and similarly if -1 e Qn. 



□ 
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We finish our discussion of tensor products with an observation for prime powers. Let 
Km denote the complete pseudograph on M vertices {i.e. an M-chque with loops): 

Lemma 2. For m ^ 3, we have G'^m = Gg K2m-z and V^im = Fg (8) K2m-3; for p an odd 
prime and 1, we have Gpm Gp® Kpm-i and Tpm = Fp (g) Kpm-i. 

Proof. We prove the results for Fpm; the results for Gpm are similar. 

• Let n = 2"'' for m ^ 3. We have q G Qn if and only if g = 1 (mod 8). Let 
r : Z,2m ^ Zg X Z2m-3 (not a ring homomorphism) be defined by T(r) = (r', k') 
such that r — Sk' + r' for r' G {0, . . . , 7}. Then, we have a — 6 G Qn if and only if 
T(a — 6) G {1} X Z2m-3, so that r induces a homomorphism F„ = Fg (g) K2m-3. 

• Similarly, for n = for p an odd prime and m ^ 1, we have q = pk' + q' E Qn 
(for q' G {0, ... ,p — 1}, which we we identify with Zp) if and only if q' G Qp. If 
T : Z2m — > ZpXZpm-i is defined by T{q) — {q', k'), we then have a—b G Qn if and only 
if T(a — b) & Qp X Zpm-i. Thus, r induces a homomorphism F„ = Fp Kpm-i. □ 

Together with Corollary 1-1, and the fact that Kpm itself may be decomposed for any 
prime p as an m-fold tensor product Kp <S> • ■ ■ <S> Kp, the graph Gn may be decomposed 
very finely whenever n is dominated by prime-power factors p"^ for p = 1 (mod 4) . 

3 Induced paths and cycles of Gn 

Even when the graph Gn does not itself decompose as a tensor product, we may fruitfully 
describe such properties as walks in the graphs G„ in terms of correlated transitions in 
tensor-factor "subsystems" . This intuition will guide the analysis of this section in our 
characterization both of the diameters of the graphs Gn, and of the factors of n for G„ a 
perfect graph. 

As Tn is a multiplicative subgroup of Z^, we may easily show that the graphs G„ 
are arc-transitive. For any pair of edges vw,v'w' G E{Gn), the affine function f{x) — 
[w' — v'){w — v)~^{x — v)+v' is an automorphism of G„ which maps v ^ v' and w ^ w'. 
Consequently Gn is vertex-transitive as well, so that we may bound the diameter by 
bounding the distance of vertices v G V{G) from G V{G), and also restrict our attention 
to odd induced cycles (or odd holes) which include in our analysis of perfect graphs. 

Let An, Bn be the adjacency graphs of the graph Gn and the digraph r„ respectively. 
We then have An = Bn = B^ if and only if —1 is a quadratic residue modulo n, and 
An — Bn + B^ otherwise; in either case, we have An oc Bn + B^. As Bn may be 
decomposed as a Kronecker product (corresponding to the tensor decomposition of F^), 
this suggests an analysis of walks in G„ in terms of "synchronized walks" in the rings 
Zp/nj by adding or subtracting quadratic units, where one must add a quadratic unit in 
all rings simultaneously or subtract a quadratic unit in all rings simultaneously. This will 
inform the analysis of properties such as the diameters and perfectness of the graphs Gn- 
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3.1 Characterizing paths of length two for n odd 

To facilitate the analysis of this section, we will be interested in enumerating paths of 
length two in G„ between distinct vertices. Because A„ oc S„ + for all n, we have 
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B„ ™,- B^ r, 



.3=1 



( B~^ ) 
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where congruence is up to a permutation of the standard basis. Thus, we may characterize 
the paths of length two in Gn between distinct vertices r,s & in terms of the number 
of ways that we may represent s — r in the form + — and — — for some 

units a, (3 e Z^; and these we may characterize in terms of products over the number of 
representations in the special case where n is a prime power. 

Definition II. For n > and r e Z„, we let Sn{r) denote the number of solutions 
{x,y) & Qn X Qn to the equation r = x + y; similarly, Dn{r) denotes the number of 
solutions (x, y) e Qn x to the equation r = x — y. 

Thus, when —1 e and An = B^ — |(_B„ + B^), the number of paths of length 
two from to r 7^ is Sn{r); otherwise, if —1 ^ Q„, the number of such paths is 
Sn{i^) + 2Dn{r) + Sn{—r). Thus, the number of paths of length two from to r reduces 
to avaluation of the functions Sn and Dn. We may evaluate these functions for n a 
prime power, through a straightforward generalization of standard results on patterns of 
quadratic residues and non-residues to prime power moduli: 

Lemma 3. For p a prime and m ^ 1, let C^m (respectively Cp^ ) denote the number 
of consecutive pairs of quadratic units (resp. consecutive pairs of non-quadratic units) 
modulo p^, and Cp^ (respectively Cp^t) denote the number of sequences of a quadratic 
unit followed by a non-quadratic unit (resp. a non-quadratic unit followed by a quadratic 
unit) modulo p"^. For primes p = 1 (mod 4), we have 



(p-5)p" 



^P 



C: 



{p - i)p 



m— 1 



(3a) 



otherwise, if p = 3 (mod 4), we have 
{p + 



c- 



c: 



{p - 3)p"'-^ 



(3b) 



Proof. As r G Z is a quadratic residue, quadratic non-residue, and/or unit modulo 
if and only the same properties hold modulo p, the distribution of quadratic and non- 
quadratic units modulo p"* is simply that of the integers modulo p, repeated times. 



It then suffices to multiply the formulae given for C++, C+ 
Aladov [9]) by p"*-^ 



C: 



(obtained by 
□ 
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Lemma 4. Let p be an odd prime, m > 0, and r e Z^m. Ifp = 1 (mod 4), we have 

for r a quadratic unit, 
for r a non- quadratic unit, (4a) 
for r a zero divisor, 

for r a quadratic unit, 
for r a non- quadratic unit, (4b) 
for r a zero divisor, 

for r a unit, 

(4c) 

for r a zero divisor. 

Proof. We proceed by cases, according to whether r is a quadratic unit, non-quadratic 
unit, or zero modulo p: 

• Suppose r G Qn- Each consecutive pair q,q + 1 G Qpm yields a solution {x,y) — 
{r{q + l),rq) G Qpm x Qp^ to x — y = r; then we have Dpm{r-) = Cpnt. Similarly, 
each such pair yields a solution {x,y) — {rq{q + l)~^,r{q + G Qpm x Qpm to 
X -\-y — r; then Spm (r) = Cpm as well. 

• Suppose r G Z^^ \. Qpm. Each consecutive pair s, s + 1 G Z^m \ Qpm represents 
a solution in non-quadratic units to x — y — 1; these may then be used to obtain 

solutions (rx, ry) G Qpm x Qpm to rx — ry = r , so that Dpm{r) = Cpm . In the case 
that p = 1 (mod 4), the negation of a quadratic unit is also a quadratic unit; in this 
case, we have the same number of solutions {rx, —ry) G QpmxQpm to rx-\-{—ry) = r, 
so that Spm[r) = C~~ as well. 

If instead p = 3 (mod 4), we instead consider quadratic units s G Qpm such that s+1 
is a non-quadratic unit. Each such pair yields a solution {x,y) = {r{s -\- 1), —rs) G 
Qpm X Qpm to X -\- y — r; then we have a solution for each such pair s, s -|- 1, so that 

Spm (r) — Cpm . 

• Finally, suppose r is a multiple of p. The congruence x-\-y = Q (mod p) is satisfiable 
for {x, y) G Qpm x Qpm only if —x is a quadratic unit modulo p for some x G Qpm, 
i.e. ii p = 1 (mod 4). If this is the case, then every x G Qpm contributes a solution 
{x, y) = {x,r — x) G Qpm x Qpm to x + y = r; otherwise, in the case p = 3 (mod 4), 
there are no solutions. Similarly, regardless of the value of p, each quadratic unit 
X G Qpm contributes a solution {x, y) — {x,x — r) G Qpm x Qpm to x — y — r. Thus 
Dpm{r) = |(p — 1) for all p; Spm{r) — |(p — 1) for p = 1 (mod 4); and Spm{r) = 
for p = 3 (mod 4) . □ 



Spm (r) — Dpm (r) 





5)p' 


lip- 


l)p^ 


lip- 


l)p^ 



rn—l 
m—l 
m—l 



for p = 3 (mod 4), we instead have 



|(p-3)p— 1 
0, 

Dr,m[r) = < 

" ^ ^ ^ i(p-l)p-^ 
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Corollary 4-1. diam(Gpm) ^ 2 forp an odd prime and m > 0; this inequality is strict if 
and only if p = 3 (mod 4) and m — 1. 

Proof. Clearly for p = 1 (mod 4) we have diam(G'pm) — 2; suppose then that p = 3 
(mod 4). Wc may form any zero divisor s = pk as a difference of quadratic units x G Qpm 
and X — pk E Qp^, so that diam(Gpm) ^ 2. We have diam(G'pm) = 1 only if is the only 
zero divisor of Z^m ; this implies that m = 1 , in which case Tpm = , so that the converse 
also holds. □ 

In Lemma 4, n = 3™" and n = 5™ are cases for which there do not exist paths of length 
two from zero to any quadratic unit. This does not affect the diameters of the graphs Gsm 
or G^m for m > 0; however, using the following Lemma, we shall see that this deficiency 
affects the diameters of G„ for any other n a multiple of either 3 or 5. 

Lemma 5. For n > odd and r e we have Sn{r) = if and only if at least one of 
the following conditions hold: 

(z) n is a multiple of 3, and r ^ 2 (mod 3); 
{a) n is a multiple of 5, and r = ±1 (mod 5); or 
{Hi) n has a prime factor pj = 3 (mod 4) such that r e 

Similarly, we have Dn{r) = if and only if at least one of the following conditions hold: 

{i) n is a multiple of 3, and r ^ (mod 3); or 
(a) n is a multiple of 5, and r = ±1 (mod 5). 

Proof. For r e Z„ arbitrary, let (ri, r2, . . . , r*) = c(r). By the decompositions = 

-Bp^mj ® ■■ ■ Bp^mt and BnB^ = Bp^m^Bj^r^^ (g) ■ • ■ Bp^mtBj^mt, we may express Sn(r) 
and Dn{r) as products over the prime-power factors of n, 

t t 

Sn{r) = n ' = n^pr^(^^) ■ 

3=1 j=l 

These are zero if and only if there exist 1 ^ j ^ t such that Spj^j {rj) = or Dpjnj [vj) = 0, 
respectively. By Lemma 4, Sp^^^j (rj) = if and only if either rj is a zero divisor of Zp^»"j 
for a prime factor pj = 3 (mod 4), or if pj G {3,5} and rj is a quadratic unit modulo 
Zp^.mj ; similarly, DpP^j (r^) = if and only if pj = 3 and rj is a unit modulo 3, or pj = 5 
and rj is a quadratic unit modulo 5. □ 

3.2 Diameter of Gn for odd n 

For odd integers n, characterizing the diameters of Gn involves accounting for "problem- 
atic" prime factors of n (those described in Lemma 5), which present obstacles to the 
construction of short paths between distinct vertices: 
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Theorem 6. Let n > 1 odd. Let 73(71) — 1 if n is a multiple of 3, and 73 (n) = 

otherwise; d^^n) = 1 if n has prime factors pj = 3 (mod 4) for pj > 3, and 6s{n) = 
otherwise; and 75(71) = 1 if n is a multiple of 5, and 75(71) = otherwise. Then, we have 

1, if n is prime and n = 3 (mod 4); 

2, if n is prime and n = 1 (mod 4); 
2, if uj{n) = 1 and n is composite; 

^ 2 + 73(71)^3(71) +75(71), ifuj{n)>l. 

In particular, diam(G'„) ^ 4. 

Proof. The diameters for a; (71) = 1 are characterized by Corollary 4-1: we thus restrict 
ourselves to the case u!{n) > 1. 

We have diam(G„) ^ 2 if and only if either Sn{r), S'„(— r), or Dn{r) is positive for all 
r & Zn ^ Tn- By Lemma 5, Dn{r) > for all r G Z„ if n is relatively prime to 15; then 
diam(G„) = 2, and r = u — u' for some u, u' G Qn for any r G Z„ if 73 = 75 = 0. If n 
is a multiple of 5, however, we have Sn{r) = Sn{—r) = D^ir) = for any non-quadratic 
unit r = ±1 (mod 5), of which there is at least one (as n is not a power of 5): thus 
diam(G„) ^ 3 if 75(72) = 1. 

Suppose that n is relatively prime to 5, and is a multiple of 3. Again by Lemma 5, 
there are walks of length two from to r if r = (mod 3), as we have Dn{r) > in 
this case. However, if n has prime factors pj > 3 such that pj = 3 (mod 4), there exist 
r G Pj'Zn such that r ^ (mod 3), in which case we have 5'„(r) = 5'„(— r) = -D„(r) = 0. 
Thus, if 73(n) = S^^n) = 1, we have diam(G„) ^ 3. Otherwise, if ^3(71) = 0, we have 
either S'„(r) > in the case that r = 2 (mod 3), or S'„(— r) > in the case that r = 1 
(mod 3). In this case, every vertex r 7^ is reachable by a path of length two, so that 
diam(G„) = 2 if 73(71) = 1 and Ss{n) = 75(71) = 0. 

Finally, suppose that either 75(71) = 1 or 73(71) = (^3 (71) = 1: from the analysis above, 
we have diam(G„) ^ 3. For r G let (ri, . . . ,r„) = cr(r), where we arbitrarily label 
P3 = 3 if n is a multiple of 3, and = 5 if n is a multiple of 5. We may then classify the 
distance of r G V{Gn) away from zero, as follows. 

• Suppose that n is a multiple of 3 and some other pj = 3 (mod 4), and that either 
n is relatively prime to 5 or r ^ ±1 (mod 5). By Lemma 5, we have Dn{r) > 
if r = (mod 3), in which case it is at a distance of two from 0. Otherwise, for 
r = ±1 (mod 3), let s = r + u for u G Qn- then s = (mod 3). Then -Dn(s) > 0, 
in which case r — u" — u' ±u for some choice of units u', u" G Qm so that r can be 
reached from by a walk of length three. 

• Suppose that 7i is a multiple of 5 and that r ^ (mod 5). We may select coefficients 
Uj G Qpf>-3 such that r^ — G {2, 3}, and such that Uj 7^ for any pj ^ 7. Let 
u = a~^{ui, . . . , Ut): by construction, we then have r—u = ±2 (mod 5) and r—u ^ 
(mod Pj) for pj ^ 7. Then either Sn{r — u) > 0, Sn{u — r) > 0, or Dn{r — u) > 
(according to whether or not 7i is a multiple of 3, and which residue r has modulo 
3 if so) : r can then be reached from by a path of length three. 
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• Suppose that n is a multiple of 5, and that r = (mod 5). If n is not a multiple of 

3, or if r = (mod 3), then D„(r) > 0; r can then be reached from by a walk of 
length two. We may then suppose that n is a multiple of 3 and r = ±1 (mod 3). 
If we also have r ^ (mod pj) for any pj = 3 (mod 4), one of Sn{r) or 5'„(— r) is 
non-zero; again, r is at a distance of two from 0. Otherwise, we have r = (mod pj) 
for any pj = 3 (mod 4), so that 5'„(r) = 5'„(— r) = -D„(r) = 0; then r has a distance 
at least three from 0. As well, any neighbor s = r±u (for u G Qn arbitrary) satsifies 
s = ±1 (mod 5). Then each neighbor of r is then at distance three from in 
from which it follows that r is at a distance of four from 0. 

Thus, there exist vertices at distance four from if ^3{n)63{n) +^5(11) = 2; and apart from 
these vertices, or in the case that 'y3{n)63{n) + 75(77,) = 1, each vertex is at a distance of at 
most three from 0. Then diam(G„) = 2 + 73(77)^3(77) + 75(77) iiu!{n) > 1, as required. □ 

3.3 Restricted reachability results for n coprime to 6 

We may prove some stronger results on the reachability of vertices from in Gn for 77 
odd: this will facilitate the analysis of perfectness results and the diameters for 77 even. 

Definition III. For a (di-)graph G, the uniform diameter udiam(G') is the minimum 
integer d such that, for any two vertices v,w & ^{G), there exists a (directed) walk of 
length d from v to w in G. 

Our interest in "uniform" diameters is due to the fact that if every vertex v G V^(r„) can 
be reached from by a path of exactly 0? in F^, then v can also be reached from by a 
path of any length i ^ d as well, which will prove useful for describing walks in F„ to 
arbitrary vertices in terms of simultaneous walks in the digraphs I'pf-j ■ 

We may easily show that F„ has no uniform diameter when 77 is a multiple of 3. For 
any adjacent vertices v and w such that w — v E Qn, we have w — v = 1 (mod 3) by that 
very fact. Then, there is a walk of length i from v to w only if £ = 1 (mod 3); similarly, 
there is a walk of length i from w to v only if £ = 2 (mod 3). For similar reasons, F„ has 
no uniform diameter for 77 even. However, for 77 relatively prime to 6, r„ has a uniform 
diameter which may be easily characterized: 

Theorem 7. Let n = p"''^ ■ ■ -p^^ be relatively prime to 6. Then 



Proof. We begin by characterizing udiam(r„), where 77 = p"* for p ^ 5 prime, using 
Lemma 4 throughout to characterize 5'„(r) for r G Z„. 

• If p = 1 (mod 4) and p > 5, we have Spm{r) > for all r G Z„; then udiam(r„) = 2. 




2 , if n is coprime to 5 and Vj : pj = 1 (mod 4); 

3 , if n is coprime to 5 and 3j : pj = 3 (mod 4); 

4 , if n is a multiple of 5. 
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• If p = 3 (mod 4) and p > 5, we have S'pm(r) = if and only if r G Z„ is a zero 
divisor. In particular, udiam(r„) ^ 3. Conversely, as |Zpm| > p™^^, there exists 
z G Qpm such that r — 2; is a unit; then there are quadratic units G Q^m, such 
that r — z = a; + I/, so that udiam(r„) = 3. 

• If p = 5, we have u G Q^-m if and only if -u = ±1 (mod 5); then r can be expressed as 
a sum of k quadratic units r = -Ui + ■ ■ ■ + -Ufe if and only if r can be expressed modulo 
5 as a sum or difference of k ones; that is, if r G {— /c, —k + 2, . . . , A; — 2, /c} + 5Z5m 
(which exhausts Zsm for /c ^ 4). 

For n not a prime power, we decompose r„ = V^^x ® ■ ■ ■ ® ^pi^t ; then a vertex r = 
cr~^(ri, . . . , rj) is reachable by a walk of length t in r„ if and only if each G ^(rp^^i) 
are reachable by such a walk in their respective digraphs. Thus, the uniform diameter of 
the tensor product is the maximum of the uniform diameters of each factor. □ 

The uniform diameter r„ happens also to provide an upper bound on distances between 
vertices in G^, under the constraint that we may only traverse walks wqW\ . . . W£ where 
the "type" of each transition wj -> Wj+i is fixed to be either a quadratic unit or the 
negation of a quadratic unit, independently for each j. More precisely: 

Lemma 8. Let n = p^^ ■ ■ -p^* be relatively prime to 6, and £ ^ udiam(r„). For any 

sequence Si,...,S£ G {0,1}, these exists a sequence of quadratic units Ui,...,U£ G Qn 
such that r = {-lY'^Ui + {-iy^U2 H h {-lY'ue. 

Proof. We first show that there are solutions to r = ui — U2 ± U3 ± ■ ■ ■ ± U£, where all 
but the first two signs may be arbitrary. We prove the result for £ — udiam(r„); one may 
extend to £ > udiam(r„) by induction. 

• Suppose n is coprime to 5: then for any r G we have Dn{r) > 0, so that there 
exist u, u' G Qn such that r = u — u' . In the case that n also has prime factors pj = 3 
(mod 4), consider s = r for any u G Qn- as there are solutions to s = u — u' for 
u, u' G Qm there are also solutions to r = u — u' ± u" . 

• Suppose n = 5"*ip^^ • • -pf*. 

— If r ^ ±1 (mod 5). Let s G Z„ be such that s = (mod 5), and s ^ 
(mod Pj) for any pj ^ 7. Then r — s ^ ±1 (mod 5), so that Dn{r) > 0; by 
Lemma 5, there are then quadratic units Ui, U2 G Qn such that r — s = Ui —U2- 
We also have 5'„(s), 5'n,(— s), /^^(s) > by construction, which can be used to 
obtain decompositions s — ±^3 ± for ^3,^4 G Qn depending on the choices 
of signs; we then have r — u^ — U2^u^^ W4. 

— If r = ±1 (mod 5), consider (ri, . . . , r^) = c(r). We select coefficients Uj, u'j G 
Qp^nj as follows. We set u[ — —ui — ri, so that 

(ri — 2mi) = (r2 + 2^2) = {ri — ui + u'l) = ±3 (mod 5). (6a) 
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For each pj ^ 7, we require Uj ^ 2 ^Vj and u'j ^ {—2 ^rj, Uj — rj}, but may 
otherwise leave Uj unconstrained; we then have 

(r-j - 2uj), {r-j + 2u'j), {rj - Uj + u'j) ^ for pj ^ 7. (6b) 

Let u = (7~^{ui, . . . ,Ut) and u' = a~^{u'i, . . . ,u[). By construction, we then 
have Dn{r — 2u), Dn{r + 2m'), Dn{r — u + u') > hy Lemma 5. There then 
exist u", u'" G Qn such that 

r —u — u +u-\- u , or (7aj 
r ^u'" -u" + u-u' = u'" -u" -u' + or (7b) 
r ^u'" -u" -u' -u' , (7c) 

selecting u" , u'" according to the desired signs for the latter two terms. 

Thus, there are solutions to r = -Ui — -^2 i ± ■ ■ ■ ±-u^ for Uj G Q„ and ^ = udiam(r„), for 
arbitrary choices of signs and r G It follows that we may decompose r = ±mi±- • ■±U£ 
for arbitrary choices of sign, provided not all signs are the same. By considering walks 
in r„ of length udiam(r„) from to either r or — r, we also have decompositions r = 
Ui + • — \-Ui and r = —Ui — • • ■ — Uj^ ior suitable choices of z^i, . . . , G Qn- □ 

The principal motivation for Lemma 8 is to bound the diameters of graphs G„ over tensor 
decompositions of the ring Z„: 

Lemma 9. Let M,N 1 be relatively prime integers, and let n — MN . Then we have 
diam(G'n) ^ max{diam(G'Ar), diam(GM)}- Furthermore, if M is coprime to 6, we have 
diam(G„) ^ max{diam(Giv), udiam(rM) + 1} as well. 

Proof. Let p : Z„ — > © '^m be the natural isomorphism. Let r G Z„ be arbitrary, 

and (r', r") = p(r). If r = {-ly^^ui H h (-1)*^^^ for some £> and Ui,...,uee Qn, 

we also have 



r' = {-iy'u[ + ■■■ + 
r" = {-lY'u'i + ••• + 



I 1 



■1) 



(8a) 
(8b) 



where {u'^u'-) — p{uj). For £ — diam(G'„), it follows that £ ^ diam(GA/) and £ ^ 
diam(GAr). 

Suppose further that M is relatively prime to 6: then udiam(rM) is well-defined by 
Lemma 8. For any a G Zjv, let £ > be the length of a walk in from to i: there are 

then ui,. . . ,ue G Qn and si, . . . , G {0, 1} such that a = (— l)''^^ H h (— 1)*^m^. If 

£ ^ udiam(rM), then for any b G Zm, there also exist quadratic units u'(, ■ ■ ■ ,u" G Qm 
such that b = (— l)*^-u'/' +■ ■ ■+ (— l)^'^^^" . We may always obtain such a walk of length i ^ 
udiam(rM) in Gjv by taking the shortest walk from to a in Gjv, and repeatedly adding 
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closed walks of length two to the end until we obtain a walk of length £ ^ udiam(rM)- 
For such a walk, we then have 

i=i i=i 

for some choice of quadratic units Uj = p~^{u'pu'-) G Qn and Sj G {0, 1}. If diam(Gjv) > 
udiam(rM), this construction yields path- lengths udiam(rM) ^ ^ ^ diam(Gjv); if instead 
udiam(rM) ^ diani(G'iv), we obtain paths of length at most udiam(rM) + l, which is satu- 
rated if there exist vertices a G V{Gn) whose distance da from is such that udiam(rM) — 
da is odd. In either case, we have diam(G„) ^ max{diam(Gjv), udiam(rM) + 1}- □ 



3.4 Diameter of Gn for n even 

The notable differences between the cases of n odd and n even are due to the sparsity of 
the quadratic units in 2,2^. compared to that of powers of other primes, and also that the 
sum or difference of two units (quadratic or otherwise) is necessarily a zero divisor if n 
is even. This results in a significant increase of the maximum diameter in the case of n 
even, compared to n odd: 

Theorem 10. Let n > even. Let 5^{n) = 1 if n has prime factors pj = 3 (mod 4) for 
Pj > 3; and S^ln) — otherwise. Then we have 



diam(G„) 



12, 


ifn 


is 


a multiple of2A; 


6, 


ifn 


is 


an odd multiple of 12; 


5, 


ifn 


is 


a multiple of 10, but not of 12; 


4, 


ifn 




8K for K > coprime to 15; 


+ 53(n). 


ifn 




6K for K > coprime to 10; 


+ S3{n), 


ifn 




AK for K > 1 coprime to 30; 


3, 


ifn 




2K for K > 1 coprime to 30; 


2, 


ifn 




4; 


1, 


ifn 




2. 



(10) 



In particular, with Theorem 6, we have diam(G„) ^ 12 for all n, and diam(Gj; 
any prime p and m ^ 0. 



^ 4 for 



Proof. We use Lemma 9 to reduce the task of characterizing diam(G'„) for n even to 
a small collection of representative cases, by factoring n = NM for suitable choices of 
coprime factors N and M. 
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• Suppose n is a multiple of 12. We may let M be the largest factor of n which is 
coprime to 12, and N — n/M. 

— li N ^ T'^r' for m ^ 3, we then have e Qat if and only \lu=\ (mod 8) 
and u = \ (mod 3), or equivalently if m = 1 (mod 24). Then Tjv consists of 
those q G 1^ such that r = ±1 (mod 24). The distance of a vertex in Gjy 
from is then characterized by its residue modulo 24, in which case we may 
show that diam(GAr) = 12. 

— Otherwise, = 4 • 3™', in which case u ^ \i and only \iu= \ (mod 4) and 
u= \ (mod 3), or equivalently \iu= \ (mod 12). Then Tjv consists of those 
5 e Zjv such that r = ±1 (mod 12); similarly as in the case above, we then 
have diam(Gjv) = 6. 

Because diam(G'M), udiam(rM) ^ 4, we then have diam(G'n) = diam(G'Ar) by 
Lemma 9. Thus diam(G„) = 12 if A?" is a multiple of 24; otherwise we have 
diam(G„) = 6. 

• Suppose n is a multiple of 10, but not of 12: specifically, n is not a multiple of 60. 

Let M be the largest factor of n which is coprime to 30, and iV = n/M . We may 
show that T]si contains only residues which are equivalent to ±1 modulo 10: 

— If n is an odd multiple of 30, we have N = 2 ■ ■ 5™'. Then w e Qtv if and 
only if u is odd, u = 1 (mod 3), and u = ±1 (mod 5); equivalently, ii u = 1 
(mod 30) or -u = 19 = —11 (mod 30). 

— If n is a multiple of 10 but not of 30, then without loss of generality = 2"^! 5™^. 
We may show that r E Qn if and only if both r = ±1 (mod 5), and 

{1 (mod 2) if mi = 1; 
1 (mod 4) if mi = 2; 
1 (mod 8) if mi ^ 3. 

In each case, we have u e Qn ii and only if w e {1, 9} (mod N) for N — 10, 
TV = 20, or TV = 40 respectively. 

As is a multiple of 10 in either case, vertices r e Zjv such that r = 5 (mod 10) can 
only be reached by a path from with length at least five, so that diam(GAr) ^ 5. 
We may show that this bound is tight by showing that every even residue can be 
formed as a sum of four elements of T„. Let x =^ y denote equivalence of two 
integers (or sets of integers) modulo m. Then, we may easily verify that 

{±1±1±1±1} =30 {26,28,0,2,4}, 
{±1 ± 1 ± 1 ± 11} =30 {8, 10, 12, 14, 16, 18, 20, 22} , 

11a 

-11 - 11 - 1 - 1 =30 6, ^ ' 

11 + 11 + 1 + 1 =30 24, 
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which proves the claim for n an odd multiple of 30. For n not a multiple of 30, 

Tjv is the set of elements q & 'L]^ such that g = 5 ± 4 (mod N) or q = —5 ± 4 
(mod N). It then suffices to show that all residues modulo 40 are exhausted by 
sums or differences of four such residues: we have 

{ (5 ± 4) + (5 ± 4) + (5 ± 4) + (5 ± 4) } {4, 12, 20, 28, 36} , 
{ (5 ± 4) + (5 ± 4) + (5 ± 4) - (5 ± 4) } {34, 2, 10, 18, 26} , 
{ (5±4) + (5±4)-(5±4)-(5±4) } {24,32,0,8,16}, 

{ (5 ± 4) - (5 ± 4) - (5 ± 4) - (5 ± 4) } =40 {14, 22, 30, 38, 6} . 

As every odd residue modulo N is adjacent to an even residue, it follows that every 
vertex in can be reached by a path of length at most five; then diam(Gjv) = 5. 
As M is coprime to both 3 and 5, we have diam(GM) = 2 and udiam(rM) ^ 3; thus 
diam(G'n) = 5 by Lemma 9. 

• Suppose that n = 8K for K coprime to 15. Let M be the largest odd factor of 
N, and = n/M = 2^ for /c ^ 3. By construction, M is coprime to 6, so that 
udiam(rM) ^ 3. We have u G Qn H and only if m = 1 (mod 8): as every odd 
residue modulo 8 can be expressed as a sum of three terms ±1, and every even 
residue modulo 8 can be expressed as a sum of four terms ±1 (with 4 requiring 
at least this many), it follows that diam(Gjv) — 4. By Lemma 9, it follows that 
diam(G„) = 4 as well. 

• In the remaining cases, we either have n = 2K for K coprime to 15 and not a 
multiple of 4, or n = 6K for K coprime to 10. We trivially have diam(G'„) = | for 
n G {2,4}; otherwise, n has odd zero divisors. As all walks of length one from in 
Gn end at quadratic units, and all walks of length two from end at even elements 
of Z„, we require walks of length at least three from to reach odd zero divisors in 
Z„. Thus, diam(G„) ^ 3. 

Let M be the largest factor of n which is coprime to 30. By construction, M is 
coprime to 6, so that udiam(rM) — 2 + S3{M) —2 + Ss{ri). 

— li n — 2K ior K coprime to 15 and not a multiple of 4, M is simply the largest 
odd factor of n, in which case n = 2*^ for A; G (1, 2}. We then have N G (2, 4}, 
so that diam(Gjv) ^jN ^2. 

— If n = 6K for K coprime to 10, we have A^ = 2 ■ 3^ for some k ^ 1. Then 
ti G Qat if and only if -u = 1 (mod 3) and is odd; that is, if m = 1 (mod 6). 
In particular, T/v contains only elements which are equivalent to ±1 (mod 6); 
from this we may easily show diam(GAr) = 3. 

In either case, it follows by Lemma 9 that 

3 ^ diam(G„) ^ udiam(rM) + 1 = 3 + ^3(71) 
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If 63(11) = 0, we then have diam(G„) = 3; we may then restrict our attention to the 
case 63(11) — 1. 

Let p : Z„ — > Zm © ^at be the natural isomorphism. As M is coprime to 15, we 
have Dm{o) > for every a G Ijm by Lemma 5: then we may express any a e 
as a difference a = m'^ — u'2 for u'^iu'2 G Qm- 

— If A?" = 2, any even residue r may be expressed as r = p~^(a, 0) = p~^(-u'i, 1) — 
p~^{u2, 1), which is a difference of the two quadratic units Uj = p~^{u'j, 1). 
Thus every even residue can be reached in G„ by a path of length two from 0. 
As every odd residue is adjacent to an even residue, we may reach any vertex 
by a path of length at most three; then diam(G„) = 3. 

— If A" = 4 or A^ is a multiple of 6, we have u G Qn if and only ii u = 1 
(mod A"), where A" = 4 if A^ = 4, and A" = 6 otherwise. We may easily show 
that the only residues r ^ Zn which may be expressed as a difference of two 
quadratic units are those such that r = (mod A"); and for any residue a = 
(mod pj), we have 5'M(a) = Sm{—(i) = by Lemma 5. Therefore, no residue 
r = p~^(a, ±2) G Z„ can be reached by a path of length two from in G„. As 
any sum of the form ±ui±U2±Us will be odd for ui,U2,U3 G Qn, such residues 
r are in fact at a distance at least four from zero. As diam(G„) ^ 3+53(71) = 4, 
it follows that diam(G„) = i = 3 + 6s{n) in this case. 

In each case, the diameters agree with the formula in (10). □ 
3.5 Perfect ness 

A graph G is perfect [10] if, for every induced subgraph H C G, the size uj{H) of the 
maximum clique in H is equal to the chromatic number x{H)- We may easily identify 
two classes of quadratic unitary graphs which are perfect: 

Lemma 11. For n even or n — p^ for p = 3 (mod 4) prime, Gn is perfect. 

Proof. If n is even, Gn is bipartite, in which case cu{Gn) = x(G„) = 2 for any non-empty 
subgraph of G„ . Otherwise, suppose that n — p^ for a prime p = 3 (mod 4) . Two vertices 
are adjacent in Gpm if and only if their residues modulo p differ. For any H C G^m, a 
maximum-size clique in H is then any set of vertices having different residues modulo p, 
where the number of different residues represented is chosen to be maximal; at the same 
time, any minimum colouring of G^m maps each residue class modulo p to a common 
colour, with different residue classes having different colours. Thus, oj{H) — x{H) for all 
such H, so that Gpm is perfect for p = 3 (mod 4). □ 

A perfect graph G contains no odd holes (induced cycles of length 2A; -|- 1 for /c > 1, 
which have no cliques larger than two but are not bipartite). Chudnovsky, Robertson, 
Seymour, and Thomas [11] characterized perfect graphs in terms of odd holes, proving a 
conjecture of Berge [12]: 
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Strong Perfect Graph Theorem. A graph G is perfect if and only if neither G nor its 
complement G contain odd holes. 

There exist imperfect quadratic unitary Cayley graphs G„ . In particular, quadratic 
unitary Cayley graphs Gp for p = 1 (mod 4) are also circulant Paley graphs, which 
Maistrclli and Penman [13] show are imperfect by exploiting the fact that they are self- 
complementarity.^ The Strong Perfect Graph Theorem implies that these graphs contain 
odd holes (again by self-complementarity) . We may consider how to "hft" odd holes in 
such graphs Gp to obtain odd holes in graphs Gn having such factors p; and we may use 
a similar strategy for Gn having distinct prime factors pi,p2 = 3 (mod 4). That is: 

Theorem 12. Gn is perfect if and only if n is even, or n = j/" forp = 3 (mod 4) prime. 

Proof. Building on Lemma 11, we demonstrate that Gn has an odd hole if n is odd and 
is not a power of a prime p = 3 (mod 4). We proceed by constructing a simpler graph 
G^, for u a factor of n, which has an odd hole. Let n — pT^^p™^ ■ • -pT* '■ 

1. Suppose that n has a prime factor which is equivalent to 1 (mod 4). Without loss of 
generality, pi = 1 (mod 4); we then let u = pi. By [13], G,y is then not perfect, and 
so Gi, or its complement has an odd hole. Note that Gj, is self-complementary, as 
multiplication of any pair of adjacent vertices by a non-quadratic unit r yields a non- 
adjacent pair, and vice-versa. Thus, G^, contains an odd hole Uq Ui U2 • ■ • u^-i Uq 
of some odd length £. 

2. Suppose instead that n has two distinct prime factors equivalent to 3 (mod 4). 
Without loss of generality, pi,p2 = 3 (mod 4); we then define v = piP2 , and let 
p : — > © be the natural isomorphism. We demonstrate the existence of 
an induced five-cycle Uq Ui U2 U3 Uq in G^, by constructing appropriate structures 
in the digraphs F^^ (g) Tp^ — p(r^) . Let 

uo = =p-\0,Q), 

We proceed by cases: 

• Suppose 2 is not a quadratic residue modulo either pi or p2 . Without loss of 
generality, we may suppose 3 ^ pi < P2; in particular, p2 — 11 or p2 ^ 19. 
By [14], there are then consecutive triples q,q + l,q + 2 G Qp.^. Let q G Qp.^ be 
a minimal such residue: as 2 ^ Qp,^ , it follows that g 7^ 1, so that g — 1 G —Qp2- 
In particular, there is an arc from g -|- 1 to 2. Define the vertices 

p{u2) = (2, g + 1) , pim) = (0, 2) , p{u^) = (2, -q) . (14) 

The vertices are distinct; in particular, U2 differs from M4, because our choice 
that 1 G implies that q + 1 ^ —q. By agreement of arc-directions Xj 

^This is in fact the simplest case of a more comprehensive theorem, which shows that the only perfect 
Paley graph is that on nine vertices. 



m = 1 =p~\i,i). (13) 
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Xk e -E'(rpi) and i/j — >■ yjt G E(rp^) for pairs of coefficients {xk,yk) — p{uk), we 
have the arc-relations 

Uq ^ Ui ^ U2 ^ U3 U4 -T- Uq (15) 

in Fj, . There are no other arc relations among the vertices Uk ■ In particular, uq 
and U3 are non- adjacent as they have the same residue modulo pi , and similarly 
for U2 and M4 ; and there are no arcs between the other pairs of vertices Uj and 
Uk because we have e ±Qpi if and only if yj - yk e tQp2 ■ 

• Suppose that 2 is a quadratic residue modulo only one of pi or p2 . Without 
loss of generality, we may suppose 2 e Qp^ , in which case P2 ^ 7 . Define the 
vertices 

p(m2)=(2,2), pM=(0,3), p(m4)=(1,4). (16) 

Clearly the vertices Uj are distinct. By agreement of arc-directions Xj — )■ G 
E(Tp-^) and yj yk & -^(^^2) for pairs of coefficients {xk,yk) = pi^k), we have 
the arc-relations 

lio — >■ lii — >■ 1*2 — >■ — >■ ii4 •<— lio (17) 

in r,^ . There are no other arc relations among the vertices Uk ■ In particular, uq 
and Us are non-adjacent as they have the same residue modulo pi , and similarly 
for ui and U4 ; and there are no arcs between the other pairs of vertices Uj and 
Uk because we have Xj Xk e ±Qpi if and only if yj - yk & l^Qp2 ■ 

• Suppose that 2 is a quadratic residue modulo both pi and p2 ■ In particular, 
Pi,P2 ^ 7, so that there arc pairs of consecutive quadratic residues modulo 
Pi and P2 (see Lemma 3). Let qi — l,qi G Qpi and 52 — 1,Q'2 G Qp^ be the 
largest such pairs in each case; in particular it follows that q2 7^ —1, so that 
—q2 — l e Qp2 . Define the vertices 

P(«2) =(^1,-^2), pius) =(0,1), p{u4) =(1,^2). (18) 

The vertices are again distinct; and by arc agreement xj — )■ Xfc G £'(rpj and 
Vj Uk ^ -^(^^2) pairs of coefficients {xk,yk) — pi^k), we have the arc- 
relations 

1^0 — >■ lil — >■ 1*2 •<— 1*3 — >■ ti4 •<— Uq (19) 

in Ti, . There are no other arc relations among the vertices Uk ■ In particular, 
Uq and U2 are non-adjacent as they have the same residue modulo pi , and 
similarly for ui and U4 ; the vertices ui and U3 have the same residue modulo 
P2 ■ By construction, we have 2^2 € Qp2 and 1 — gi G —Qp^ , so there is no arc 
between U2 and U4 ; and similarly there is no arc between uq and U2 ■ 
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In each case, G^, contains an induced cycle Uq Ui U2 ■ ■ ■ U(_-i Uq for some odd £ ^ 5. Let 

A* = { (20) 

' I m^ mo 'c ^ ^ 

we may then obtain a similar odd hole in by identifying each Uj with a corresponding 
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G {0, . . . , — 1} C Z^. Then — Xk & if and only if — Uk & Qu for any 
1 ^ j, k ^ i, which implies that xq xi ■ ■ ■ xe-i xq is a cycle without chords in G^. 

Let A?" be the largest factor of n which is coprime to /j, {i.e. N = n/ii), and let 
T : Z„ — > © Ztv be the natural isomorphism. 

• Suppose u = pi = 1 (mod 4). If n is a multiple of 5, we may suppose that pi = 5 
without loss of generality; then N is coprime to 5. By Lemma 5, we then have 
SN{r) > for any — r e Q^: then Gn contains a closed walk r for some 
u E Qn- By concatenating this walk with |(£ — 3) copies of a closed walk of length 
two at 0, we obtain a closed walk yo Z/i ■ • ■ Vi-i Ho in Gjy . We may then construct 
a walk 

C=(xo,yo) {xi,yi) ■■■ (x^_i,y^_i) (xo,yo) (21) 
in Gfj, (g) G]\f: because Xq Xi ■ ■ ■ X£_i Xq is an induced cycle, so is C. 

• Otherwise, suppose fi = p^^p^^- If n is a multiple of 3, we may suppose that 
Pi = 3 without loss of generality; then N is coprime to 3. By Theorem 7, we then 
have i > udiam(rAr), in which case by Lemma 8 we may construct a closed walk 
yo yi ■■ ■ ye-i yo in Gn by setting yo = 0, and letting yj+i - yj e ±Qn whenever 
Xj+i — Xj G {i.e. with agreement in the signs), and similarly for yo — Vi-i- 
Define the walk C in Z^ © Zjy as given in (21): we then have {xj^i^yj^i) — {xj, yj) G 
±(Q/. ® Qi^) 6ach j, and similarly {xo,yo) - {x£_i,yt-i) G ±((5^ © Qat). 

In either case, if we define vertices Vj = T~^(xj,yj) G V{Gn), the walk Vq Vi ■ ■ ■ Vf_i Vq is 
an induced cycle of odd length in Thus Gn is not perfect, unless n is even or a power 
of a prime p = 3 (mod 4). □ 

4 Decomposing symplectic operators mod n 

Our final result is a bound on the complexity of decompositions of symplectic operastors 
modulo n, which follows from the bound on the diameter of We may define the 
symplectic form (modulo n) as the 2m x 2m matrix 



(22) 



the symplectic group modulo n '&Y>2m{'^n) is the set of 2m x 2m linear operators S {sym- 
plectic operators) with coefficients in Z„ such that S~^a2mS — (72m- 
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Convention. For operators S e Sp2mi'^n) for a fixed m, we will adopt the convention of 

indexing the rows and columns by integers modulo 2m, starting with 1. Thus, for a row 
k G {m + 1, . . . , 2m} in the "bottom" half of a matrix S, the row A; + m G {1, . . . , m} will 
be in the "top" half, and vice-versa. 

Symplectic operators are clearly invertible operations, and therefore may be reduced 
to by Gaussian elimination. We also consider a variant procedure, in which row- 
operations are constrained to themselves be symplectic. For the operator definitions 
below, actions of operators are defined via the action of left-multiplication on a square 
matrix over Z„. 

Definition IV. For row- indices j, e {1, . . . , 2m}, a symplectic row operation acting on 
Sp2^(Zn) is one of the operators Mj"'\ Ej^k , Cj^k , or C~l defined as follows: 

• For any a e Z^, let /xj""* e GL2n(Z„) be the hnear operator which multiplies the j*^ 
row of its operand by a. Then, we define Mj"^ — /j,^"^ l^%m ■ 

• Let Ej^k G GL2m(Zn) be the linear operator which exchanges rows j and k of its 
operand. Then we define 

p / iij-k = m (mod 2m); 

Sj^k£j+m,k+m otherwise, 

for A*j~^^ as defined above. 

• Let Xj,k G GL2m(Zn) be the linear operator which adds row j of its operand to row 
k. Then we define 

{Xj,k , \i i -k = m (mod 2m); 

X3,k Xllm,3+m ' if 1 ^ J, ^ ^ Or m + 1 ^ j, /c ^ 2m; (24) 
Xj,k Xk+m,j+m , Otherwise. 

These operations are defined so as to be symplectic themselves; we wish to demonstrate 
an upper bound to the number of such symplectic row operations required to transform 
an arbitrary symplectic operator to the identity. 

Hostens et al. [15] provide a decomposition of symplectic operators into 0(m^ log(n)) 
symplectic row operations, in an application to the the decomposition of an important 
family of unitary operators for quantum computation (specifically, the Clifford group over 
qudits of dimension n). We refine this decomposition to obtain an upper bound to O(m^), 
giving an upper bound which is independent of the modulus n. 
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4.1 Reduction to greatest common divisors modulo n 

We first describe the decomposition of [15] in detail. Tlie main concept is to reduce 
S e Sp2^(Z„) to anotfier operator 5" wfiicfi acts trivially on, e.g., the standard basis 
vectors e^, e2m- This reduces the problem to decomposing an operator S e Sp2^_2(^n), 
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(25) 



Embedding the matrix groups Sp2(Z,„) C ■■■ C Sp2m_2(^n) ^ ^P2mi'^n) in the manner 
described above, one may recursively apply this process to obtain a sequence of symplectic 
row operations which multiply to transform S to /2m- As the inverse of each symplectic 
row operation is also a symplectic row operation, this yields a decomposition of S. 

The reduction from S to S' as above is performed by a hybrid of Gaussian elimination 
and Euclid's algorithm for computing greatest common divisors. We illustrate this on 
a 2m X 2 matrix [ v w ] , for a pair of column vectors v = [ vi f 2 • ■ • V2m ] ^ and 
w — [wi W2 • • • W2Tn V subject to the constraint w^(j2mV = 1. By performing suitable 
symplectic row- additions, we may simulate the Euclidean algorithm in the second column, 
for each pair of rows {j, j + m) for j e {1, . . . , m}, to obtain 
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(26a) 



computing "greatest common divisors" (modulo n) in the lower block in the second col- 
umn, and using these to clear the upper block. We then perform further row-additions to 
compute further greatest common divisors in the second column, in pairs of rows {j, j + 1) 
for j e {m -|- 1, . . . , 2m — 1}, in to perform the following transformation of the the second 
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column: 
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(26b) 



Note that as w^(72mV = 1, there is an integer combination of the coefficients of w which 
is equivalent to 1 modulo n; then 

gcd(7i, . . . ,7„) = gcd{wi,...,W2rn,n) = 1. (27) 

The above row-transformations then transform the two-column matrix [ v w ] as follows: 



Vl 
Vm. 
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Vm+2 Wm+2 
V2m W2m 
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(28) 



where v'^ = 1 follows from e^„jO'2„tv' = "w~^o'2m^ = 1- We may repeat the sequence of 
transformations to compute greatest common divisors in the first column, for each pair of 
rows {j,j + 171) for J e {1, . . . , m}, and subsequently in row-pairs + 1) in the upper 
block: 



gcd((/?i,(/?2) 
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(29) 
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The reduction of [15] applies this procedure for v = S'e^, w = Se2m- Applying these 
transformations to S yields a matrix S' as illustrated in (25), as the other columns S'sk 
for k ^ {m, 2m} must satisfy 

^~kS'a2m^m = hSa2mSem = 0, (30) 

and similarly ^ S' (J2rn^2rn = 0. The complexity of a single iteration of this reduction 
is 0(m log(n)), which arises from the cost of repeating Euclid's algorithm (expressed in 
fixed-width integer addition steps) 0{m) times to reduce the m*'^ and 2m^^ columns to 
Bm and e2m respectively. Iterated m times over all column-pairs 6^,6^+^, we obtain the 
upper bound of 0(m^log(n)) reported by [15]. 



4.2 Improved upper bounds via the diameter of Gn 

The complexity of the above decomposition may be reduced to O(m^), by substituting 
an explicit simulation of Euclid's algorithm via symplectic row transformations with a 
product of constant size. This is possible by using short paths in the graphs Gn to reduce 
the number of addition steps in order to obtain coefficients 'jj and Lpj (or coefficients 
equivalent to them, up to a multiplicative unit) using a constant number of row-operations. 

The primary obstacle to reducing the complexity of a single iteration of the reduction 
of [15] is the computation of greatest common divisors in row-pairs + m), arising 
from constraints on obtaining "derived" row-additions on these row-pairs. The iterated 
operator (for a G {0, 1, ... ,n — 1}, which we identify with a G Z„) can be easily 
obtained in constant depth ioi j ^ k + m (mod 2m) and a G , by the equality 

Cr,k = Mf'^ C,,, , (31a) 

which one may verify by the action on standard basis vectors. However, C^"^-^^ cannot 
be decomposed in this manner: the closest we may come is in the case where a — for 
some ti G , in which case we have 

Cr,,-C^l^ Mf^ ) . (31b) 

We may apply the result of Theorem 10 as follows: 

Lemma 13. For distinct row-indices j, /c G {1, . . . , 2m} and for any a G there exists 
a sequence of units Oi, . . . , G and signs Si, . . . , G {—1, +1} for some i ^ 12, such 
that 

= Wfc^l Mf"' Cll Mf"' ■ ■ ■ Mf'^'-'Cli . (32) 

Proof. It suffices to note that as diam(G'„) ^ 12, there exists such a sequence of signs 
and quadratic units Ui, . . . ,Ui G Qn such that a = SiUi + S2U2 + • • • + siUi. We may 
then take either aj = Uj (in the case that k ^ j + m) or a unit aj such that Uj = a| (in 
the case that k = j + m), and apply the decompositions of (31) to obtain the desired 
decomposition.^ □ 



■^For k ^ j + m, we may in fact obtain the further bound of £ ^ 3, as the diameter of the unitary 
Cayley graph X„ = Cay(Z„,Z^) is at most three [4]. 
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We may apply this to reduce the complexity of decomposing symplectic operators as 
follows. We use the following additional Lemma, whose proof is deferred to the appendix: 

Lemma 14. Let 7 = gcd(x, y, n): then there exist a,b,c & Z such that ax + by + cn — ^ 
and where both a and b are relatively prime to n. 

For a vector x = [ X2 ■■■ X2m Y i let 7^ = gcd(a;j, a;j+,„, n) for each j e 
{1, . . . , m}. Let aj be coefficients such that aj+rnXj+m + cijXj = 7j (mod n) as guaranteed 
by Lemma 14, and define rj — aj^^aj: we then have 
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(33) 



Each coefficient 7^ generates the same additive subgroup as 7^ modulo n; ii di, . . . ,dm are 
coefficients such that Xj — dmOj^^rnlj — dmlj-, we then have 



C, 



m+1,1 



C2mTm [Xl ■■■ Xm 71 • • • 7m ] ^ = [ ■ ■ ■ 7i • • ■ 7™ ] ^. 



(34) 



The above performs the reduction of (26a), up to multiplicative units, in 0(m) symplectic 
row operations. We may similarly emulate the reductions of (26b) and (29) in Oim) 
symplectic row operations, using Lemma 14 to reduce the computation of greatest common 
divisors (up to multiplicative unit factors) to performing powers of the operators (7,,^. 

To summarize, using the bound on the diameter of the quadratic unitary graph G„, we 
may refine the decomposition of symplectic operators in [15] by substituting an explicit 
simulation of Euclid's algorithm by a constant-size sequence of symplectic operations. 
This substitution provides an upper bound of O(m^) for a decomposition of an operator 
S e Sp2m(Z„), a bound independent of the modulus n. 



5 Remarks and open problems 

It should be noted that quadratic unitary graphs, while easy to describe, arc closely tied to 
unsolved problems in computational complexity theory. In particular, testing adjacency 
in a graph Gn is precisely the quadratic residuacity problem, which has no known efficient 
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algorithms and is considered unlikely to be efficiently solvable (see e.g. Chapter 3 of [16]).^ 
Because of this, an efficient algorithm (deterministic or randomized) for discovering the 
shortest path between two vertices in Gn should be considered unlikely. We may then 
ask whether there are efficient algorithms for discovering "short" paths (having length 
bounded by a fixed constant) between vertices in Gn- 

In Section 3.5, we provided a partially non-constructive proof that odd holes arise in 
quadratic unitary graphs G„ which are odd but not a power of a prime p = 3 (mod 4). 
Numerical investigation suggests that, in particular, five-holes (odd holes of size five) are 
very common in those Gn which are not perfect graphs, even when restricting to five-holes 
involving the arc — > 1. It would be interesting to obtain a classification of all five-holes 
which occur in the imperfect graphs Gn ■ 

As we noted in the introduction and in Section 3.5, the graphs G„ for n = 1 (mod 4) 
prime are also Paley graphs. Shparlinski [18] shows that prime-order Paley graphs these 
graphs have high energy {i.e. the operator 1-norm of the adjacency matrix), coming to 
within a factor of (1 — ^) of the upper bound Euia,x{n) = ^n{y/n + 1) shown in [19] for 
graphs on n vertices. We may ask to what extent this and other properties of circulant 
Paley graphs generalize for quadratic unitary graphs. 

The author would like to thank Charles Matthews, Robin Chapman, and Chris Godsil 
for helpful discussions. This work was partly written during a visit to School of Computer 
Science at Reykjavik University. This work was performed with financial support by 
MINOS EURONET and the EURYI scheme. 
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A The existence of special Bezout coefficients 



For a sequence of integers Xi,X2, ■ ■ ■ ,Xk, Bezout coefficients are a corresponding sequence 
of integer coefficients Oi, 02, . . . , Ofc sucfi tfiat gcd(a;i, . . . , Xk) = ^ the existence of 
sucfi a sequence of coefficients Oi, . . . , is implied by the "simple" Euclidean algorithm. 

Consider the greatest common divisior of a sequence xi, . . . ,Xk together with another 
integer n: this is equivalent to computing 7 = gcd(a;i, . . . , x^) modulo n via Euclid's 
algorithm. We may compute greatest common divisors modulo n recursively, by comput- 
ing 72 = gcd(a;i,a;2) modulo n, then 73 = gcd(gcd(a;i, 0:2), 0:3) modulo n, and so forth. 
However, for each intermediate stage 1 < j < A;, it is not necessary to obtain 7^ itself, but 
instead a similar residue 7^ which generates the same subgroup modulo n; by definition, 
the set of integer combinations modulo n of such an integer is the same as the set of 
integer combinations of jj, so that gcd(a, 7^) = gcd(a,7j) (mod n) for any a G Z. 

The simplest application of this observation is that in any integer x may serve as 
a substitute for its own greatest common divisor with n: 

Lemma 15. Let 7 = gcd(a;, n) for x,n e Z; then there exist a,b & such that ax+bn — ^ 
and where a is relatively prime to n. 

Proof. Consider arbitrary a,b & such that ax + bn — 7. Let a — gcd(a,7): then a 
divides a, x, and n. We have -x = (- — b—) G -Z. By construction, - is an integer 
relatively prime to - : thus, a; is a multiple of -. Let m = — : then — = - G Z, and 
furthermore is relatively prime to n. Then, if we let a G Z be such that ^ = 1 (mod n), 
we have ax = 7 (mod n) as required. □ 

Wc generalize the above lemma as follows. In order to compute a suitable integer 7^ 
which generates the same additive group (modulo n) as 7^ = gcd(7j_i, Xj) for each j, we 
may compute Bezout coefficients a, 6, c such that 



where aa = l (mod n) , in which case we may let 7^ = 67^. That is, if such a G Z^ exists, 
we may compute 7^ as the sum of Xj with some multiple of 7j-i, which can be computed 
using a single addition operation and a single scalar multiplication. We show that such 
Bezout coefficients may always be found by proving Lemma 14 (page 23): 

Lemma 16. Let 7 = gcd(a;, n): then there exist a,b,c & Z such that ax + by + cn = j 
and where both a and b are relatively prime to n. 

Proof. Let x' = gcd{x,n) and y' = gcd{y,n): by Lemma 15, we then have x = Uxx' 
(mod n) and y = Uyy' (mod n) for multiplicative units Ux,Uy G Z^, and 7 = gcd{x',y'). 
Define 



axj + b^j-i + cn = 7^. 
If we may find such a set of coefficients that a is coprime to n, we then have 

Xj + alffj^i = a7j (mod n) , 



(35) 
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these are both divisors of n, and form a relatively prime pair. We may then partition 

the prime factors of n into those which divide x, those which divide y, and those which 
divide neither. Let be the largest factor of n whose prime factors divide x, Ny be the 
largest factor of n whose prime factors divide y, and Nn = n/N^Ny : then A^^; and Ny are 
coprime, so that Nn is also an integer and relatively prime both to N^ and Ny. We then 
have n — N^NyNn- 

As X and y are coprime, there exist integers a,b E 1^ such that ax + by = 1. Note that 
a is coprime to y , from which it follows that a is coprime to Ny as well, as Ny and y have 
the same prime factors; similarly, b is coprime to N^ . Let 

if gcd(a, N^Nd) = gcd(6, NyNn) = 1; 
if gcd(a, N^Nd) > 1, but gcd(6, NyN^) = 1; 
if gcd(a, 7V,7V^) = 1, but gcd(6, 7V^7V„) > 1; ^ ' 

otherwise; 

and let a — a + hy and (3 — b — hx. If a has prime factors in common with N^Nn, then a 
does not, by the fact that both y and yNy are relatively prime to N^Nn ; otherwise, a is 
coprime to N^Nn anyway by the coprimality of a to N^Nn. In either case, we also have a 
coprime to Ny, by the coprimality of a and y. Thus, a is relatively prime to n = N^NyNn] 
and similarly, /3 is coprime to n. We may then observe that 

ax -\- j3y — (a + hy)x + (6 — hx)y — ax + by — 1 , (39) 

from which it follows that ax' + f3y' — 7. Let a, 6 e Z be such that u^d = a (mod n) and 
Uyb = P (mod n) : then, we have 

ax + by = ax' + (3y' = 7 (mod n) : (40) 

as a, Ux, (3, and Uy are all coprime to n, both a and 6 are also coprime to n. □ 
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